Breaking of Chiral Symmetry and Spontaneous Rotation in a Spinor Bose-Einstein 

Condensate 
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We show that a spin-1 Bose-Einstein condensate with ferromagnetic interactions spontaneously 
generates a topological spin texture, in which the m = =bl components of the magnetic sublevels 
form vortices with opposite circulations. This phenomenon originates from an interplay between 
ferromagnetic interactions and spin conservation. 



The Mermin-Ho texture jl| in superfluid 3 He describes 
an interesting thermodynamic equilibrium state, in which 
a circulation remains nonvanishing in a cylindrically sym- 
metric vessel, despite the vessel being at rest. This phe- 
nomenon is due to the fact that the boundary condition 
at the surface of the vessel imposes a topological structure 
on the I vector. In this Letter, we show that the spon- 
taneous formation of such a topological spin texture also 
occurs in a Bose-Einstein condensate (BEC) of atomic 
gases with spin degrees of freedom 0, . 

In contrast to the Mermin-Ho texture in 3 He, the phys- 
ical origin of the spin texture formation proposed in this 
Letter is the interplay between ferromagnetic interactions 
and spin conservation. Consider a spin-1 BEC with fer- 
romagnetic interactions in an m = magnetic sublevel. 
The spin-exchange collisions between the atoms transfer 
the m = state into m = ±1 states, as + — > 1 + (— 1). 
As a consequence, magnetization in the x-y plane may 
arise due to the ferromagnetic nature of the interaction. 
However, uniform magnetization of the entire system is 
prohibited because of spin conservation, which results in 
various spin textures [3, la, H ? such as staggered domain 
structures. In this Letter, we show that the topological 
spin texture is spontaneously generated as a result of the 
spin exchange dynamics under spin conservation. 

In this topological spin texture, each of the m = ±1 
components contains a vortex, whose directions are oppo- 
site. Therefore, from the symmetry of the Hamiltonian, 
there are two degenerate textures: one has + and — vor- 
ticies in the m = 1 and —1 components, respectively, 
and the other has — and + vortices. We show that the 
symmetry between these two textures is spontaneously 
broken in the course of the dynamics even when the ini- 
tial state has chiral symmetry. 

Topological spin structures in a spinor BEC of an 
atomic gas have been realized by the MIT group |g, 
who obtained topological phases adiabatically imprinted 
on the spin components by a quadrupole magnetic field, 
producing coreless vortices. The stability of such topo- 
logical spin structures has been studied by several au- 
thors 0, El EH • Recently, it was predicted [H 03 
that the dipolar interaction creates a coreless vortex state 
through a mechanism similar to the Einstein-de Haas ef- 



fect. 

We consider a system of spin-1 Bose atoms with mass 
M confined in a potential V. The Hamiltonian of the 
system is given by |3, |4| 
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where ijj m is the field operator for an atom in the mag- 
netic sublevel m = 0, ±1, H = -ft 2 V 2 / \2M) + V, 
and F is the spin-1 matrix. The spin-independent and 
spin-dependent interactions are characterized by go = 
47rh 2 (a + 2a 2 )/(3M) and g 1 = 47rft 2 (a 2 - a )/(3M), re- 
spectively, where as is the s-wave scattering length for 
the scattering channel with total spin S. In the case of 
spin-1 87 Rb, ch is negative and the ground state is ferro- 
magnetic 0,E3- 

When the potential V is axisymmetric with respect to 
the z axis, the Hamiltonian is invariant under spatial 
reflection with respect to an arbitrary plane containing 
the z axis, e.g., (x,y) — > {—x,y). This transformation 
changes a clockwise vortex oc e~ 1 ^ into a counterclock- 
wise vortex ex e 1 ^ with azimuthal angle (j). Hence, if only 
one of them is spontaneously realized, we call it chiral 
symmetry breaking. Also, from the symmetry of the 
Hamiltonian, the x, and z components of the total 
spin and the z component of the orbital angular momen- 
tum are conserved independently. 

We consider the case in which the initial state is in the 
m = mean-field ground state satisfying 



(#o+#o|^o| 2 )^o = Mo*o. 



(2) 



If the m = ±1 components are exactly zero, is a 
stationary state of the multicomponent Gross-Pitaevskii 
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(GP) equations, 
ih^- = (Ho +g n)il>o + ^=(F + il> 1 +F-il>- 1 ), (3a) 
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function, n 



where ip m is 

the m = ±1 components is analyzed by the Bogoliubov 
de Gennes equations, given by 



= |^i | 2 - |V-i| 2 , and F+ = F*_ 
i). The stability against excitation in 



[J?o-Mo + (5o+ffi)|*o| 2 ]4V5i*o4t = e^u% 

(4a) 

[//o-Mo + (ffo+5i)l*o| 2 ]^t+ffi*u 2 ^ = -e ( Ht, 

(4b) 

where and are the eigenfunctions for a Bogoli- 
ubov mode with eigenenergy From the axisymmetry 
of the system, the Bogoliubov modes can be classified ac- 
cording to the angular momentum i, for which u 



oc e" 



and Vm oc e 



\ We find that u^p couples with v^\* in 



Eq. (gj, and hence the excitation of the m = 1 component 
with vorticity £ is accompanied by the m = — 1 compo- 
nent with vorticity — £, as a consequence of the orbital 
angular momentum conservation. If all the eigenenergies 
are real, the state is dynamically stable. If there ex- 
ist complex eigenenergies, the corresponding modes grow 
exponentially and the state \I>o is dynamically unstable. 

For simplicity, we restrict ourselves to two-dimensional 
(2D) space. This situation can be realized by a tight 
pancake-shaped potential V = Mu 2 (x 2 + y 2 + X 2 z 2 )/2 
with A ^> 1, where the axial confinement energy is so 
large that the dynamics in the z direction are frozen. 
In this case, the interaction strengths can be charac- 
terized by the dimensionless parameters g 2B = 
gj[X/(27r)] 1 ^ 2 N/(ay lo huj), where N is the number of 
atoms, aho = [hf(Muj)] 1 ^ 2 , and j = 1, 2. 

We numerically solve Eq. (0) by the imaginary-time 
propagation method and diagonalize Eq. Q to obtain the 
Bogoliubov spectrum for the state \I/o- Figure shows 
the lowest Bogoliubov energies for £ — and ±1 as a 
function of g 2D , where is determined by #o D A£? D = 
go/ gi — —216.1 which is the ratio for spin-1 87 Rb [l(j(- In 
the parameter regime shown in Fig. d the three modes 
exhibit complex eigenenergies. A crucial observation is 
that there is a region (—3.9 > g 2B > —10.7) in which only 
e^ 1 ) are imaginary. This indicates that these two modes 
only are dynamically unstable in this region, where one 
mode has vortices oc e ±l( ^ and the other mode has vortices 
ex e^ %< ^ in the m = ±1 components. These two modes are 
degenerate because of the chiral symmetry of the system. 
In this region, we expect that the rn = ±1 components 
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FIG. 1: (Color) Real and imaginary parts of the lowest Bo- 
goliubov energies for £ = 0, =bl, where the label £ indicates 
that the m = ±1 components of the eigenfunction are pro- 
portional to e ±ie< ^. The two energies e^ 1 ^ are degenerate due 
to the axisymmetry of the system. We have taken the param- 
eters of spin-1 87 Rb atoms, where the spin-independent inter- 
action strength is related to the spin-dependent strength 
fl P by fl0 2D = -216.1 5l 2D 



start to rotate, despite there being no external rotating 
drive is applied to the system. 

In order to confirm the spontaneous rotation phe- 
nomenon predicted above, we numerically solve the GP 
equation © in 2D using the Crank-Nicholson scheme. 
The interaction strengths are taken to be g^P = 2200 and 
g 2B — —10.18, with the ratio gQ B /g 2B again chosen to be 
that of spin-1 87 Rb. For this set of interaction parame- 
ters, Ime^ /(hu) = 0.0707 and all the other Bogoliubov 
energies are real. The initial state is the ground state of 
Eq. (|2J) for the m = component plus a small amount of 
random noise in the m = — 1 component. To extract the 
Bogoliubov excitations from -0 m (t), we define ^3 
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which represents the degree of excitation in the modes 
up 1 ^ oc v^ 1 " 1 * oc e ±l( ^. Figure [21 shows the time evolution 
of P±i and the density-phase profile of each component 
at cut = 130. We find that P±i increase according to 
exp^Ime^^/ft), with their initial ratio P-i/Pi ~ 7.5 
kept constant, resulting in an exponential growth in the 
angular momenta of the m = ±1 components. Thus, the 
m = ±1 components spontaneously rotate if the initial 
noise has angular- momentum fluctuations. 

Suppose that we can prepare an initial state of the sys- 
tem in which the chiral symmetry is preserved to great 
accuracy, say, P-\jP\ = 1.0002. Then, the result in 
Fig. indicates that vortices will not be created in the 
time scale in which the linear stability analysis is ap- 
plicable. For a longer time scale, however, the chiral 
symmetry is spontaneously broken due to the nonlinear 
effect. Figure 01 (a) shows the time evolution of the frac- 
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FIG. 2: (Color) Degree of Bogoliubov excitation P±i in 
Eq. 0. The interaction strengths are g 2 ^ — 2200 and 
g 2r> = —10.18. The initial state is the ground state of 
Eq. (0> for the m — component plus a small amount of 
random noise in the m — — 1 component. The dashed line is 
proportional to e °- 141ujt . The insets show the density-phase 
profiles at out = 130, where the size of the frame is 16 x 16 in 
units of aho = [Ti/ (Mou)] 1 ^ 2 . The heights of the vertical lines 
in the insets indicate \iprn\ 2 CLho/N, which is 0.001 for m = ±1 
and 0.0025 for m = 0. 

tion n_i = J dr\ip-i\ 2 /N and the orbital angular mo- 
mentum per particle L \ = —i f drip*L -^d^-i/ {Nri-i) 
for the m — —1 component. The m = 1 compo- 
nent is given by n\ ~ n_i and Li ~ — L_i from spin 
and orbital angular momentum conservation. The ini- 
tial value for the m = — 1 component is taken to be 
= 10- 4 r(e^ + 1.0001e-^)^ , which gives P-i/Px ~ 
1.0002. As long as this ratio is kept constant, the forma- 
tion of the vortex states shown in the insets of Fig. 01 is 
not expected. In fact, as shown in Fig. 01(b), no vortex is 
created around the first peak of n_i at out ~ 100. How- 
ever, at out ~ 160, L-i starts to deviate from (blue solid 
curve) and the chiral symmetry is dynamically broken. 
Consequently, the vortex states emerge in the m — ±1 
components at out c± 340 as shown in Fig. 01(c). 

The instability in the state with chiral symmetry 
[Fig. 01 (b)] against forming the rotating state [Fig. 01 
(c)] implies that the energy of the latter is lower than 
that of the former. In order to confirm this, we take the 
energy dissipation into account by replacing i on the left- 
hand side of E q. (Ot with i — 7 [l8| . The time evolution 
with 7 = 0.03 PJj is shown by the red curves in Fig. 01 
(a) and clearly indicates that the energy of the rotating 
state [Fig. 01 (d)] is lower than that of the state having 
chiral symmetry [Fig. 01(b)]. For 7 = 0, L_i oscillates 
with a large amplitude due to the excess energy released 
from the initial state, while for 7 = 0.03, the sign of the 
angular momentum is unchanged. 

The bottom panels in Figs. 01 (b)-(d) show the spin 
vector distributions. In Fig. 01(b), the magnetic domains 
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FIG. 3: (Color) (a) Time evolution of the fraction n_i 
(dashed curves) and the orbital angular momentum per par- 
ticle L-i (solid curves) in the m = — 1 component with 
(7 = 0.03, red) and without (7 = 0, blue) dissipation. The 
interaction strengths are the same as in Fig. The initial 
state is given by = ^ 0, ^-1 = 10" 4 r(e^ + l.OOOle _i0 )tf o , 
and ipi — 0, where ^0 is the ground state solution of Eq. J3J. 
(b)-(d) Snapshots of the density-phase profiles and spin tex- 
tures. The heights of the vertical lines in the density-phase 
profiles show \ipm\ 2 Q>ho/N, which is 0.001 for m = ±1 and 
0.0025 for m — 0. The length of the vector in the bottom 
panels is proportional to (F 2 + F 2 ) 1 ^ 2 and the color repre- 
sents (F 2 + F y 2 ) 1/2 /(|^-i| 2 + |^o| 2 + |^i| 2 )- The size of the 
frame is 16 x 16 for the density-phase profiles and 12 x 12 for 
the spin textures in units of aw 



in the opposite spin directions are separated by a domain 
wall at x = 0. On the other hand, topological spin struc- 
tures are formed in Figs. 01 (c) and (d). The underlying 
physics of the spin structure formation is the interplay 
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FIG. 4: (Color) Energy and density-phase profiles obtained 
by the variational method with Eq. JSJ). The interaction 
strengths are the same as in Fig. The size of the frame 
of the insets is 16 x 16 in units of aho and the heights of the 
vertical lines show \ipm\ 2 a^ /N, which is 0.0005 for m = zbl 
and 0.0025 for m = 0. 



between the ferromagnetic interaction and spin conser- 
vation. The growth in the spin vectors must be accom- 
panied by spatial spin structure formation to conserve 
the total spin angular momentum. It should be noted 
that the area in which the length of the spin vector is 
long is larger in Figs. |3| (c) and (d) than in Fig. |3| (b), 
since the spin vectors must vanish at the domain wall in 
the latter. This is why the energy of the state in Fig. [3| 
(d) is lower than that of Fig. [31(b). That is, the formation 
of the topological spin structure increases the (negative) 
ferromagnetic energy of the system more than does the 
formation of the domain structure. 

The above energy argument concerning the spin do- 
main and topological structures can be reinforced by ap- 
plying the variational method. We assume the variational 
wave function to be 






where ^ 



V 







1^1 



(6) 



l^-il 2 ) 1 / 2 so that the total 



density n is kept to be |^o| 2 irrespective of the values 
of the variational parameters, reflecting the fact that the 
spin-exchange process hardly changes the total density 
because go ^> \g±\ for spin-1 87 Rb atoms. We minimize 
the energy of the system calculated from Eq. ((SJ with 
respect to c, x, an d x' f° r a given as shown in Fig. 21 
The state at = 7r/4, which is similar to the state shown 
in Fig. |3|(b), has a maximum energy, while the topolog- 
ical spin states at — and ir/2 have minimum energy, 
in agreement with the above discussion. 

The results presented above can be realized using cur- 



rent experimental setups. For example, when the radial 
trapping frequency is uj = 100 x 2tt Hz and its ratio to the 
axial trapping frequency is uj/uj z = 0.01, the interaction 
parameters for Figs. 12141 correspond to N ~ 8880 spin-1 
87 Rb atoms. The time scale for the appearance of the 
topological spin structure (e.g., out ~ 300 for the initial 
condition in Fig. |3J) is ~ 0.5 s. If the ratio \gi/go\ can 
be increased with a decrease in go by Feshbach resonance 
or by using other atomic species, we can reduce the time 
scale, e.g., to about 1/10 for \gi/go\ — 1- 

In the presence of an external magnetic field the 
linear and quadratic Zeeman terms enter the Hamilto- 
nian (^J. Since the total spin is conserved, the linear 
Zeeman term only rotates the spin at the Larmor fre- 
quency and does not affect the dynamics. When the 
magnetic field is applied in the z direction, the quadratic 
Zeeman effect raises the energy of the m = ±1 com- 
ponents compared with the m = component. If the 
quadratic Zeeman energy exceeds the ferromagnetic en- 
ergy, the m = state becomes the ground state and no 
excitation to the m = ±1 components occurs, which is 
the case for B > 400 mG for the parameters uj = 100 x 2tt 
Hz, uj/uj z = 0.01, and N ~ 8880. We have numerically 
confirmed that the dynamics in Figs. and |3| are quali- 
tatively unchanged for a magnetic field of ~ 100 mG. 

In conclusion, we have proposed a novel mechanism 
of spontaneous formation of a topological spin structure 
in the spin-1 BEC prepared in the m = state. The 
m = ±1 components increase exponentially from initial 
random seeds due to dynamical instabilities and form 
singly-quantized vortex states (Fig.[2J). Even if the clock- 
wise and counterclockwise rotation components are as- 
sumed to be equal in an initial seed, one of them even- 
tually becomes dominant (Fig.[3J). This chiral symmetry 
breaking is attributed to the fact that the topological spin 
structure is energetically the most favorable due to the 
ferromagnetic interaction. This spontaneous spin struc- 
ture formation is essentially caused by the spin exchange 
dynamics under the constraint of spin conservation. We 
expect that many more interesting spin textures may also 
be spontaneously generated in isolated spinor BECs. 
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